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Abstract 

We flow a hypersurface in Euclidean space by mean curvature flow 
with a Neumann boundary condition, where the boundary manifold is 
any torus of revolution. If we impose the conditions that the initial 
manifold is compatible and does not contain the rotational vector field 
in its tangent space, then mean curvature flow exists for all time and 
converges to a flat cross-section as t —> oo. 
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1 Introduction 

We consider Mean Curvature Flow (MCF) of hypersurfaces with a Neumann 
boundary condition, choosing the boundary manifold to be an n-dimensional 
torus of rotation of any profile embedded in (n + l)-dimensional Euclidean 
space. If the initial manifold is compatible with the boundary condition 
and is transversal to the rotational vector field, then the flow exists for 
all time and converges to a flat cross section. A tool regularly used with 
MCF with a Neumann boundary condition is to assume convexity of the 
boundary manifold [1] [9], which gives a sign on boundary derivatives and 
allows application of Hopf maximum principle. Clearly we cannot impose 
this in the case of the torus. Instead we observe that we are essentially 
considering a graphical problem and use a Stampacchia iteration argument 
similar to that used by Huisken [5] for graphical mean curvature flow. 

Mean curvature flow with a Neumann boundary condition has been con- 
sidered as graphs in the perpendicular case by Huisken [5], and also for 
more general angles in dimension 2 by Altschuler and Wu [lj. In both cases 
the flow exists for all time and converges to a special solution: In the first 
case to a flat plane, while in the second to a translating solution. The level 
set method in the case of the right angled Neumann condition in a convex 
cylinder has been studied by Giga, Ohnuma and Sato [12] . Considering 
the perpendicular angle condition further, Stahl|9j[10j showed that if the 
boundary manifold is a totally umbillic surface and the initial manifold is 
convex, then under mean curvature flow the manifold shrinks to a point. 
Furthermore, renormalising homothetically about this point the flow con- 
verges to a half sphere. Buckland [2] used a monotonicity argument (again 



with a perpendicular boundary condition) to classify the Type I boundary 
singularities of a mean convex initial surface. The perpendicular Neumann 
boundary condition has also been considered in flat Minkowski space by the 
author [7J , where the boundary manifold was chosen to be a convex timelike 
cone and the rescaled flow converges to a hyperbolic hyperplane. 

Suppose S C M n+1 is a smooth orientable hypersurface with outward 
pointing unit normal vector [i. Following Stahl [10J we say 
F : M n x [0, T] — > M n+1 satisfies Mean Curvature Flow with a Neumann free 
boundary condition S if 

'% = H = Hv V(x, t) G M n x [0, T] 
< F(-,0) = M 

'F(i,t)cSl V(x, t) E dM n x [0, T] 

(i/,/z)(a?,t) = V(x, t) G dM n x [0, T] 

where v(x, t) is the normal to F at time t. We will often write M% = F(-, t). 
In this paper we choose S to be a rotationally symmetric torus of any profile 
- topologically S 1 x S 1 ™ -1 - and we flow a disk contained within the interior 
of this torus. At a point p = (pi, . . . ,p n +i) G M n+1 we define 



r = \]pl+P 2 n +n r = -(°! • • • ,0,p n ,p n+1 ), G = -(0,...,Q,-p n+ i,p n ) . 

We prove the following: 

Theorem 1.1. Suppose J] is a torus of rotation and Mq is an embedded 
disk satisfying the boundary condition which nowhere contains the vector 
field 6 in its tangent space, then a solution to equation |Ip with initial data 
Mq exists for all time and converges uniformly to a flat cross-section of the 
torus. 

The proof uses an integral iteration technique to obtain the cruicial gra- 
dient estimates. In |5j, Huisken uses similar arguments in the case of a cylin- 
drical boundary manifold. The advantage of this method is that boundary 
curvature is less of an issue; we require boundedness of the derivatives of 
certain functions as opposed to a sign on them. We remark that this theo- 
rem allows some unusual initial manifolds, for example the disc may wrap 
itself around the inside of the torus several times (see Figure [I]) as long as 
it is transversal to the vector field generated by the group of rotations. 

For this paper we will need various geometric quantities on various man- 
ifolds. A bar will imply quantities on R n+ , for example A, V, . . . and so 
on; no extra markings A, V, . . . will refer to geometric quantities on Mt the 
flowing surface at time t and for any other manifold Z, A z , V z , . . . etc. will 
refer to the Laplacian, covariant derivatives, ... on Z. We will define the 
volume form on Mt to be (i and define fig to be the volume form on dMt- 
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Figure 1: A jesters hat initial manifold is taken to a flat disc as t — > oo. 

2 The torus 

In this section we make some remarks about S, a torus of revolution. We 
define the half space = {( Xl ,...,x n ,Q)\xi el,i„>0} C M n+1 , where 
we will sometimes write y = xiei + . . . x n _ie n _i and x n = r. Suppose we 
have any compact domain Q, C with smooth boundary dVt parametrised 
by P : S n — > M™, then by rotating in the {e n , e n+ i}-plane we define X 
to be the torus of revolution d£l sweeps out. Since 8 is the direction of the 
rotation, at a point p 6 £ we know that 6(p) £ T p £. 

We will require the values of the second fundamental form of £ in the 
direction explicitly. We parametrise £ by 

3(x,6) = P(x) - (P(a;),e n ) e n + (P(x),e„) [cos#e n + sin#e n+ i] . 

If v p is the outward pointing unit normal to P in W 1 then //, the outward 
pointing unit normal to £ in M ra+1 is given by 

A* = v p - (u p ,e n ) e n + (y p , e n ) [cos#e n + sin(9e n+ i] . 

We may easily see that 

d 2 J I dP 



8*86 \d X f e " } -»'^,,+n,H0e (( + J j 

is perpendicular to /u. Therefore we know that the direction S = r6 is an 
eigenvector of A^(-, •). The eigenvalue may be calculating by writing 

d 2 3 

Q2Q = ~ ( P ' e «) t COS + Sin ee n+l] 

and so A^(r6, rO) = — (y p ', e n ) (P, e n ) = — r (fx, r) . 
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It will also be useful to consider the flowing manifold as a graph over 
SI in the interior of the torus, when this is possible. At every point x G f2 
we assign u, the angle through which we need to rotate x about to hit 
the manifold, so that we may parametrise the manifold inside the torus by 
F : Q R n+1 , where 

F(x) = y + r(cos(u)e n + sm(u)e n+1 ) . (2) 

We may now compute all standard geometric quantities with respect to 
u by standard methods. For example 

9ij = $ij + r 2 DiuDjU , 

-r(y + D n u(cos(u)e n + sin(u)e n+ i)) - sin(u)e n + cos(u)e n+ i 
v „ , 

V 

where we define v = det(gjj) = yT+ r 2 \Du\ 2 . Similarly we may calculate 
the equations for mean curvature flow in these coordinates for Mo a manifold 
which may be parametrised as above by uo. Such a parametrisation will 
always be possible when Mq is transversal to 0, although the range of no 
may be more than 2tt. Considering graphically, equation ([!]) is equivalent 
to 

u(x, 0) = uq(x) Vx G O 

^=gVD ijU +^(l + ^)=-tf VxGft (3) 
j-Du = \/x£dn 

where 7 is the outward unit normal to dfl and g 11 is the inverse of the 
metric in this parametrisation. We note that uniform parabolicity of the 
above is equivalent to the gradient estimate v < C < 00. We also remark 
that v - 1 1 



(v,— sin(u)e n +cos(u)e n+1 ) (v,0\ F ) ' 

Evolution equations, boundary derivatives and 
initial estimates 



Here we will obtain initial estimates on various quantities via a maximum 
principle of the following form: 

Theorem 3.1 (Weak Maximum Principle). Suppose we have a function 
f : M n x [0, T) ->• E then if f satisfies 

(jt _ /(p, t) < V(p, t) G M n x [0, T) such that V/(p) = 
(V/, n) < V(p, t) G dM n x [0, T) 

then f(x,t) < sup /(p,0) for all (x,t) G M n x [0,T). 

p6M" 
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We will repeatedly use the following easily verified relations 

Vr = r, W x r=(X,0)-G, V x e = -(X,0)-r . (4) 

r r 

Lemma 3.2. Let 9 be the angle around the torus, taken from some arbitrary 
base angle. Then 

dt r 



A)9=-(r T ,Ve) = --^(u,0)(u,r) , 



and at the boundary 



V„0 = O . 



Proof. Using cylindrical coordinates on M n+1 we see that V0 = f and from 
this we may calculate the evolution equation of 9. We see that 

d9 H 



and 



AO <■'" [ '-It — \ - - /— e) - — \ - /- ft 



2 'fjA-nltv 



Therefore 

At the boundary we have that 0(p) G r p S, and therefore 

□ 

Since the above only depends on the derivatives of 9, the same equations 
hold for manifolds that wrap themselves around the torus more than once 
(as in Figure [I]) by simply extending the range of 9 to be more than 2ir. The 
function u defined in the previous section is an example of this, and on the 
flowing manifold u will satisfy the same evolution equation as 9. From now 
on we use the extended 9. 

Corollary 3.3. The function 9(-,t) is bounded above by its initial supremum 
and below by its initial infimum. 

Proof. This follows directly from the above maximum principle. □ 
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Therefore, the disk may not twist itself around the torus any more than 
it is twisted initially. The following will be required later: 



Lemma 3.4. The function r evolves by 



dt- A)r 



while at the boundary 

V^r = (/x,r) 



Proof. Similarly to Lemma 3.2 we have % = —H (r, v) and calculate 



^ r ^ ^ dx l dx^) \ dx^} 



— f)\f) \ + / r. V » - V « ) 



^ \\ r \ dx l ' / ^' dxi / \ ^af» dxi ^ dxi / 

-J^<^, r) + 1 |6> T | 2 . 

r 



□ 



Remark 3.5. We will assume throughout that ro < r < r%, a consequence 
of the flowing manifold staying within the torus. Certainly this will be true 
for all time that we have a gradient estimate, and apriori will be true for 
t £ [0, e) for some small e > 0. 

We will need the following well known evolution equations: 

Lemma 3.6. On the interior of a manifold moving by mean curvature flow 
the following hold 



dgij 
dt 



-2Hhij (6) 



±-a)h = H\A? (7) 

Proof. See for example [U Lemma 3.2, Lemma 3.3 and Corollary 3.5] . □ 

The boundary derivative of H and relations on the second fundamental 
forms of the flowing manifold and the boundary manifold for equation ([I]) 
were calculated by Stahl [HI Proposition 2.2, 2.4], and are summarised in the 
following Lemma: 
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Lemma 3.7. At the boundary 

and for X G T p M n T p £ then 

A E (X,u) + A{X,n) = . 

We now define the gradient function, v = j^-^- Without loss of gen- 
erality we may assume that this is positive and so may define the related 
function Q = log v. 

Lemma 3.8. The evolution equations for Q and v (while they are finite) 
are 

(|_ A )g__ W ._|VO|» , (1-a)„ = -„|^-^ , 

and the boundary derivatives are 

V M Q = -A s (u, v) , V^v = -vA^{v, u) . 

Proof. We will first calculate the evolution of w = {rO, v). 
Using Q we may immediately see that 

^ = -H(V u (rO),v) + {r6,VH) = {r9,VH) . 

Writing Z = rO we have 

A (u, Z) = g* (WiiWju) - V v .j^v, Z) 

+ 2 g v (Wiu, VjZ) + gv (u, %Z - HV u z) . 

For the first of these terms, take a orthonormal basis {fi, . . . , f n } at a 
point p 6 M. We extend this to give orthogonal geodesic coordinates at p. 
We calculate that at p, 

9 tj (Vi(^i/) - V v ,u, Z) = g ij (Ufiu), Z) 

= g ij Vjhu9 lk (ffc, Z) - g ji h u g lk (h jk v, Z) 
= V z tH-(u,Z)\A\ 2 

where we used the Weingarten and Codazzi equations. Since the right hand 
side does not depend on the coordinate system this holds for all p & M, 
For the second term we have 

°" = o ,s ('■ w) " T - (£>•") rT > = ■ 
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The final term also vanishes; we may see that 
V 2 XY Z = V Y ((X, r)6- (X, 6) r) - V WxY Z 

= ^ [(X, 6)6- (X, r) r + (X, r) r - (X, 6) 6] = . 

r 

Therefore 

and the evolution equations for v and Q immediately follow. 

At the boundary using Lemma |3.7| and the fact that /ilflwe have 



V„{u,Z) = A(»,Z T )+(v,{r, fJ ,)0- (ji,e)T) = (u,6)(^r)-A j: (u,Z-wu) . 

Now from calculations in Section [2] we know that Z is an eigenvalue of 
j4 s (-, •), and using the calculation of this eigenvalue we see 

V„ (v, Z) = (v, 6) (ji, r) + - (v, 6) A X (Z, Z) + wA^{y, v) = wA^(v, v) 

r 

and we are done. □ 

Corollary 3.9. While v is bounded we have the estimate H 2 < Ch for some 
constant Ch > depending only on Mq and E (and independant ofv). 



Proof. We calculate 

' d_ 
di 



A \ H 2 v 2 = -2v 2 \X7H\ 2 - 8Hv (VH, Vv) - 6H 2 |Vt;| 2 

= -2 [v 2 \ VH\ 2 + Hv (VH, Vv}] 

-Q[H 2 \Vv\ 2 + Hv{VH,X7v)] . 

At a positive stationary point we have HVv = —vVH, and so the above 
vanishes and we may apply our maximum principle. Since ^ < r we have, 
using Remark 3.5, a uniform bound. □ 

Remark 3.10. We may attempt to get a positive lower bound on H in a 
similar way and indeed the evolution equations are amenable. In fact due 
to the boundary condition, this is not useful. Using the vector field Z as in 
Lemma 13.81 we see 

= — I {fJ>,Z) dfiQ = [ div(Z T )dfi = f — dfi 
JdM Jm Jm v 

where again we used equation Q. Since by assumption v is initially positive, 
H cannot be initially positive everywhere, and therefore (for example) weak 
mean convexity implies we must have a minimal hypersurface - not much 
left for the flow! Indeed a corollary of our main theorem is that the only such 
minimal hypersurface that satisfies our initial conditions is the flat profile. 
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4 Integral estimates 



As a prerequisite to applying the Stampaccia iteration method as in [5] we 
now give some of the required boundary estimates. In particular, we modify 
various Lemmas for graphs with boundary in [3] to manifolds with boundary. 

For a start, we will require the Michael-Simon-Sobolev inequality from 
[6]. While this holds in much more general situations, we will only require 
M to be smooth embedded n-dimensional manifolds in M n+1 . 

Lemma 4.1 (The Michael-Simon-Sobolev inequality). There exists a con- 
stant Cs > depending only on n such that for any function f G C 1 (M) 
such that f has compact support, we have 

n-1 

( / \f\^d(i) <C S I \Vf\ + \H\\f\dfi . 

\JM J JM 

We also need such an inequality not just on functions of compact closure, 
but functions that may be non-zero at the boundary dM. 

Lemma 4.2. For any compact manifold M with boundary dM and for any 
function f G C l (M) we have 



n-1 



/ / '<///) <Cs 

JM 



\Vf\ + \H\\f\djl+ [ \f\djj, a 

JdM 



' M JdM 

where the constant Cs depends only on n. 

Proof. This is as in [3l Lemma 1.1]. Let d : D — > K be the function giving 
the minimum distance along the manifold to the boundary. This is smooth 
close enough to the boundary. We define for k large enough rjk = min{l, kd}, 
and let % be a C 1 smoothing of this. We consider the sequence fk = Vkf 
for k G N. Since p,({x\f(x) ^ /»}) — > as i — > oo we have that 

n— 1 n—1 

( f \f k \^dfi) U -+( [ \f\^dil) n , / \H\\f k \djl^ [ \H\\f\dfi. 

\JM J \JM J JM JM 



We also see, 



M 



\Vf k \dfi< I |V/MA+ / \f\\Dr}kW 

JM JM 



The first term of the above may be estimated similarly to the other terms. 
For the final term we choose a special parametrisation of the collar. We 
parametrise by C : dM x [0, e) — > M. n+1 for e small enough by setting C(x, e) 
to be the point obtained by starting at x G dM and moving distance e down 
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the geodesic starting at x with direction — \i. Therefore = Vd and the 
metric induced by C has gi n = <5j n . Therefore for k large enough 



/ \f\\D Vk W< ( k\f\dfi 

JM J{x&M\d(x)<±} 



k / \f\\/det(gij{x,s))dfL d ds 
Jo JdM n v 

\f\ydet(gij(x,0))djj,d= / \f\dfrd 



>dM n v JdM" 
as A; — )• oo. □ 

Clearly we need to estimate boundary integrals, and we now give one 
way of doing so, based on [31 Lemma 1.4]. 

Lemma 4.3. For M a compact manifold with boundary, then for all f E 
W 1 '°°(M) we have 



\f\djb<CE [ |V/| + (|#| + l)|/|d/2 

JM 



IdM 

where the constant > depends only on E. 

Proof. This is essentially just divergence theorem. We now use d, the min- 
imum distance to S in IR n+1 and note that at S, Vd = —fJ>. We take a 
smooth function (f) : R — > M such that (j>'(0) = — 1 and </>(a;) = for x > R 
where R is less than the minimum focal distance of S. We define = 4>(d) 
- a smooth function on M n+1 . Then 

M> = g^V 2 4> (S) - H (vM) < Ci(l + • 



dx l J \dxi J 

for some C\ > depending on the derivatives of and so 
fdfig = / div(/V^)d/2 



dAf 



Ai 



(V/,V0) + /A^£ 



<C 2 / |V/| + /(|fT| + l)d/2 , 

for some C2 > depending on the derivatives of (f>. □ 

Corollary 4.4. For all f E C 1 (M) i/iere exisis a constant C$ depending 
on n and S suc/i £/iai 

rt-l 

/ " <C S [ |V/| + (|#| + l)|/|d/i 
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Due to the boundary condition we have the following: 

Lemma 4.5. Suppose f : M n x [0, T) — > IR is once differentiable in time 
such that 4z,f € L 1 (Mt). Then the following holds for t > and f3 = 0: 

d f /" d / rr2. 



/<# = / f-H'fdjl (8) 

Proof. The perpendicular boundary condition implies that the manifold 
does not flow out through the boundary. Therefore, we know that in the 
parametrisation defined by F in ([!]) over the stationary domain M n that 

/ fdfi= / fJdet(gij(x,t))dx . 
J Mt JM n v 

Now the equation follows from □ 



Remark 4.6. Integrating equation (4.5) for / = 1 with respect to time and 
rearranging we see that 

/ H 2 dfidt <\M \ , 

JM t 

that is we have a parabolic L? estimate on H that does not depend on the 
time interval. 

We will also require the following well known Lemma, which serves to 
streamline the iteration argument of the next chapter: 

Lemma 4.7. Suppose (j) : (&o,oo) — > M. is a non-negative non-increasing 
function such that for all h > k > ko then 

where C, a and j3 are positive constants. Then if j3 > 1 then 4>(ko + d) = 
for 

d a = C[0(fc o )] /3 - 1 2 a ^i . 
Proof. See Lemma 4.1 i)], for example. □ 

5 The gradient estimate via iteration 

Here we will give a bound on the gradient Q = log(— (i/,rO)) via integral 
estimates. We define Qk = (Q — k)+ and A(k) = {x G M\Qk > 0} and we 
aim to get suitable estimates on the quantity 

\\A{k)\\ = F I djldt , 



JO JA(k) 

ultimately showing that this is zero for all sufficiently large k. We begin 
with some LP estimates on ft. 
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Lemma 5.1. There exsits a k± > such that k > k\ and even p 

fT 







M 



Q p k djldt<C Q {p)\\A{k)\\ 



where Cq > depends on p, n, E, Mq and k\. 



Proof. Using Lemma 3.2 and writing w = (rO, v) we have that 

2 



d 
dt 



A 



= \e Xd 
<\e xe 



(u,9) {u,t) - \\V9\' 



w r e 



A 



w 



for Cf > 0. At the boundary this function has zero derivative in the fx 
direction (from L emm a 3.2). On A(k), we may estimate w < e~ k and so 

l 1 9 1 ° 

and writing \V6\ = - 



recalling Remark 



3.5 



-, we may get a positive 



lower bound on |V0| for k sufficiently large. Therefore, we may choose a 
large enough A, ko > so that for all k > ko on A(k) 

A J e Ae < -C||V^| 2 



dt 



< -3d 



holds, for C|, C| > depending on A, where we used our bounds on 8 from 
Corollary 3.3 We now agree to write C n for any bounded positive constant 
depending only on Mq, p, S, ko and n. Using Lemma 3.8 we calculate that 
on A(k) for k > ko and p > 2, 



A e X9 Q p k < Q p k - 2 [-CiQt\V6\ 2 -pe Xd Q k \VQf 



d 
dt 

- e xe p(p - l)\VQ\ 2 + C n pQ k \VQ\\Ve\] 

< Q P k ~ 2 [c n Q k - C e 2 Ql\V6\ 2 - C n Q k \VQ\ 2 - C n \VQ\ 2 

< C n Q p ~ 2 - 2C e 3 Ql - C n Q p - l \VQ\ 2 - C n Q p - 2 \VQ\ 2 , 

where we used the lower bound on | V0| and Young's in equa lity of the form 

and divergence 



2 l2 

"' + %- c repeatedly. We may now use Lemma 



ab - 2 , 2e 
theorem to see that 



4.5 



d 
dt 



M 



e xe Q p k dfi < 



d 

m \ dt 



A\e X9 Q p -H 2 e X9 Q p dil 



+ 



dM 



pC n CzQ p k l dfi d 
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Estimating as above and using Lemma 4.3 and Corollary 3.9 



d_ 
dt 



[ e xe Q p k dp 

JM 








< [ qV 

JM 


C n 


- 2ClQl - C n Q k \VQ\ 2 


< [ Ql 2 

JM 


C n 


— 


dp, . 



dp, 



Hence choosing k > ki = max{/co,supQ + 1} we ma Y integrate to get 

M 



QUpdt < C 



JA(k) 



Qu 2 dpdt 



JA(k) 



It is easily verified that the above argument also holds if p = 2. There- 
fore, by induction we may estimate for p even and k > k\ 



T 



JM 



Q p k dpdt<C Q (p)\\A(k)\\ 



□ 



In addition to the above we also need an estimate on ||A(A;)|| which 
does not depend on T. This comes about by using Remark |4.6| to estimate 
the highest order terms in a Laplacian after careful use of the boundary 
condition: 

Proposition 5.2. There exists a k2 > such that for all k > ki there exists 
a constant C depending only on Mq, £ and n such that 

\\A(k)\\<C 

Proof. We consider the bounded function k = 6 2 r 2 . We will calculate the 



time derivative of the integral of / over the manifold. From Lemmas 3.2 
and 13.41 we see 



d_ 
dt 



A K = r s 



(Vr,V0) - 2|V0| 



+ 



-2|0 T | 2 -2|r T | 2 



80r(Vr,W) . 



At the boundary we have V^k = 6 2 V^r 2 and so calculate 

Vf,Kdpg= [ div(8 2 Vr 2 )dp = [ (V6 2 ,Vr 2 ) + 6 2 Ar 2 dp 

JM JM 



/■ 

JM 



48r{V6,Vr) + 6' z -2rH (u, r) + 2|0 T | 2 + 2|r 



,Ti2 



dp 
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Therefore, by divergence theorem, 

d 



d_ 
dt 



M 



M 



dt 



A I k — H 2 Kdfi + / V uHidp, 



c)M 



[ -2r 2 \V0\ 2 - 26 2 rH (u,v) - H 2 ndjl 
Jm 



We note that (v, r) 2 + (u, 0) 2 < \v\ 2 = 1 and so (v,r) 2 < l-{v,0) 2 = \6 T \ 2 . 
Since |V#| 2 = ^ 2 , using Young's inequality we see 



d 
dt 



[ K dft< [ -2\6 T \ 2 + 2(6 2 r\H\)(\{v,r)\) - H 2 fdji 

-\e T \ 2 + c\H 2 dfi . 



M 



for some C\ > by the boundedness of r and 8. 

Now integrating with respect to time as in Remark 4.6 and using the L 2 
bound on H, 



e T \ 2 diidt < d 



JM 



JM 



H z dfldt+ I 6 z r^dfl 

M 



t=o 



for some constant C3 > depending on the bounds on 9 2 , r 2 and \Mq\ but 
not on T. On the region A(k), (is, 6) < \e~ k and so choosing ki large enough 
that (v, 6) < ^ then 



\\A(k)\\<2 



JM 



\0 T \ 2 dildt <2C 3 . 



(9) 

□ 



We now put these together to give the gradient estimate. 



Theorem 5.3 (Gradient Estimate). There exists a Cq > depending only 
on n, S and Mq such that for all time Q < Cq. 



Proof. By Lemma 3.8 we may calculate for p > 2 that 



£ 

dt 



a) Ql < pQl 1 (±-a\q- p(p - l)Q p k - 2 \VQf 



<-pQ P k ~ 2 (Qk\VQ\ 2 + (p-l)\VQ\ 2 ) . 



Using C n as in Lemma 5.1 we see using the bound on \H\ and Lemma 4.3 
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then 



d_ 
dt 



M 



Q V k dfi< I Ql 2 - p Q k \VQ\- 2 -p(p-l)\X7Q\ 2 -H 2 Ql 



M 



+ C n Q k + C n \VQ\\dji 
< f Q p k ~ 2 [C n Q 2 k + C n ] dfL-C n f P Q p k - 1 \VQ\ + + l)Q p h dfi 

JM JM 



< C 2 I Ql' 2 + Q p k djl- d 



M 



rip 



M 



where on the last line we used Corollary 4.4 with / = Q?. Choosing k > 
ks = max{ sup Q, ki, k^\ and integrating with respect to time 



sup / Q p k djd + C\ 
te[o,T] Jm 



np 

QZ^dfL 



M 



dt < C 2 



JM 



Q P k ~ 2 + Q P k dfidt . 



We now deal with the left hand side as usual (see [5] ) , and so after repeated 
use of the Holder inequality we have for even p, 



p(n + l) 



Q, n dfidt < C n 



JM 



n + l 



f T f Ql" 2 + Q P k dfidt 

JO JM 



<C 3 \\A(k)\\ 



where we used Lemma 5.1 For k 3 < h < k, the Holder inequality now 
implies 



\h-k\ p \\A(h)\\< 



JM 



Q p k dfidt<Cz\\A{k)\Y 



n + l 



We may now apply Lemma 4.7 to get that = for k = k% + D where 

DP = C 3 2 2+n ||^(fci)||^ . 
Setting p = 4, by Proposition |5.2| the Theorem is proved. □ 



6 Long time existence and convergence 

We now state and prove the main Theorem. 

Theorem 6.1. Suppose E is a torus of rotation, Mq is a manifold satisfying 
the boundary condition that nowhere contains the vector field in its tangent 
space. Then a solution to equation £/p with initial data Mq exists for all time 
and converges uniformly to a flat cross-section of the torus. 
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Proof. We take $7 to be a cross-section of the torus S and rewrite the mani- 
fold as a graph, no, over the cross-section as in section [2] so that the manifold 
may be parametrised by equation Q. At a point on the flowing manifold, 
this will be equal to the function 9. As noted in section [2] for both uniform 
parabolicity of equation ^ and a gradient estimate on u we need to bound 
the function v = j^-^ . We also note that while v is finite we may write Mt 
as a graph. 



Since v = rv, Theorem 5.3 gives the upper bound v < C, and so we 
have both uniform parabolicity and a gradient estimate. Corollary |3.3| also 
gives C° bounds on u. Therefore, by standard methods we have existence 
for all time. For example since equation Q has linear boundary conditions, 
with trivial modifications we may apply the arguments of [HI Section 8.2 and 
Chapter 12]. 

For convergence we consider integrals of the derivatives of the graph over 
0. We have J* = -Hv and so using our gradient estimate and Corollary 



) dxdt= I I H 2 v 2 dfidt<d / / B 2 dfidt<C 2 



fdu\ 

JSl\dtJ Jq JM JO JM 

where C\, Ci > are constants independant of T. 

We see that in coordinates |0 T | 2 = 1 — (v. B) 2 = ^=f^ = ^ Du ) . Therefore 
using the gradient estimate again 

f [ \Du\ 2 dxdt<C 3 [ [ ^^ vdxdt = C 3 f [ \6 T \ 2 dildt <C 4 
Jo Jn Jo Jn v Jo J m 

for constants C3 , C4 > where we used equation ^ . 
Therefore there exists a constant C > such that 

+ \Du\ 2 dxdt<C . 

Writing ua(t) = iqt Jq u { x , t)dx for the integral average of u at any time then 
by the Poincare inequality, the above is enough to ensure that u(x,t) — > 
Ufi(t) uniformly as t — > 00. Since we also have from Corollary |3.3| that 
inf u(x, t) is nondecreasing and sup u(x,t) is nonincreasing then in fact 

u(x,t) converges uniformly to a constant as t — > 00. This corresponds to 
uniform convergence of Mt to a flat cross-section of the torus. □ 
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